hep-ph/0611194 



Determining Neutrino and Supernova Parameters with a Galactic 

Supernova 

Solveig Skadhauge and Renata Zukanovich Funchal 

Instituto de Fisica, Universidade de Sao Paulo, 
C. P. 66.318, 05315-970 Sao Paulo, SP, Brazil 

Abstract 

We investigate the sensitivity of some of the proposed next-generation neutrino experiments to a 
galactic supernova. In particular, we study how well the supernova parameters (the average energies 
and luminosities) can be separated from the unknown neutrino oscillation parameters (^13 and the 
neutrino mass hierarchy). Three types of experiments, all in the 100 kilo-ton class, are compared. 
These are: a 540 kton water-Cherenkov detector, a 100 kton liquid Argon detector and a 50 kton scin- 
tillator detector. We demonstrate that practically all of these proposed detectors have the possibility 
to determine the hierarchy of the neutrino masses if the angle ^13 is sufficiently large (sin^(0i3) ~ lO"'^) 
and the hierarchy of the average energies is larger than about 20%. They can at the same time de- 
termine some of the supernova parameters well. The average energy of the z/^ and species can be 
determined within 5% uncertainty in most of the parameter space suggested by supernova simulations. 
The detection of several separable channels measuring different combinations of charged current and 
neutral current processes is crucial for determining the value of ^13 and the hierarchy. However, there 
are cases where a few of the SN parameters can be determined rather well even if only the main 
charged current detection channel is available. 

PACS numbers: 14.60.Pq,25.30.Pt,97.60.Bw 
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I. INTRODUCTION 



The fortuitous observation of a handful of neutrinos from the supernova 1987A |2| started 
the era of experimental supernova (SN) neutrino physics. It confirmed our main ideas of the 
physics of a supernova from a core collapse, although also a set of (minor) disagreements were 
found. However, the low statistics of neutrino events collected from SN1987A, partially due 
to the fact that it happened at a distance of roughly 50 kpc from Earth, was not sufficient to 
really extract much information. Evidently, with a larger detector and a closer by supernova, 
the prospects for gaining high accuracy information about both supernova and neutrino physics 
are immense. In the case of a galactic SN, one can observe of the order of 500 events per kton 
of detector material and there are several proposals for neutrino detectors in the 100 kton 
range j^. 

A type II supernova is the death of a giant star and the huge emission of light as well as 
neutrinos is an effect of the gravitational core collapse of the star. In fact, about 99% of the total 
energy is emitted in the form of neutrinos and antineutrinos in a roughly 10 second interval. 
The consecutive burning of different elements, structured in an onion shell form, ends at the 
silicon burning phase. This phase produces an iron core and as fusion of iron is impossible, 
the gravitational collapse of the core is triggered once it reaches a well-known size, related to 
the Chandrasekhar limit. The core heats up during the implosion and photo-disintegration 
of the iron atoms starts, having as a byproduct free protons and neutrons. The subsequent 
neutronization gives rise to a z/g flux from the deleptonization process e + p Uf, + n. The 
SN densities in the interior of the core are so high that even neutrinos are trapped. However, 
neutrinos, interacting only weakly, manage to escape near the surface of the core region. A 
solid feature in all SN simulations is the formation of a shock wave caused by the rebounce of 
material falling on to the core. It is believed that this shock wave will cause the explosion of 
the outer material of the star. However, it has for a long time been a pending embarrassment, 
that SN simulations are incapable of producing the SN explosion as the shock wave looses too 
much energy to the media and halts before erupting the outer parts of the star. It has been 
suggested that neutrino physics might play a crucial role in re-energizing the shock wave [4|, 
but until now a successful calculation has not been performed. After this so-called collapse 
phase, where a neutron star or, in rare cases, a black hole is formed, there is an accretion phase 
and a cooling phase. 

In this article we will study the neutrino flux from the cooling phase. It is probably the 
phase where the physics is best controlled, although many things are still unknown. The 
energy spectra will be almost thermal, as the densities in the SN are so high that in fact the 
neutrinos will be in thermal equilibrium. But what are the average energies of each neutrino 
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flavor? Again, a robust feature seems to be tlie liierarcliy of tliese average energies. Electron 
neutrinos can interact in the medium both through charged current (CC) and neutral current 
(NC) exchanges, whereas t'^ and z/^ suffer only NC interactions as there are no muons or taus in 
the SN material. As a consequence, electron neutrinos will maintain their thermal equilibrium 
to a larger radius, therefore escaping with a smaller temperature or, equivalently, lower average 
energy {Ee). Furthermore, as z/^ and v-ri both, only interact via NC, their energy spectra will 
have identical properties. We will treat them as indistinguishable, denoting them with the 
common index x. Anti-electron neutrinos also interact via CC reactions but with a smaller 
cross-section (since there are less protons than neutrons) resulting in a higher temperature 
(average energy (-Eg)) than electron neutrinos. So, a hierarchy of the form {E^ < (Eg) < (E^) 
is predicted to exist. However, the exact values of the average energies as well as the strength 
of the hierarchies, vary quite substantially in different simulations and also depend on the type 
of the progenitor star [sfHi. 

Neutrino flavor transitions have been observed in atmospheric, solar, reactor and accelerator 
neutrino experiments. The simplest and most widely accepted way to explain these transitions 
is to allow neutrinos to have masses and mixings. Although, the neutrino oscillation parameters 
have been determined to increasingly astonishing precision during the past few years 
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currently, 



+ 7.3 X 10"W < Aml^ < +9.0 x 10"W 

0.25 < sin^Ou < 0.37 (1) 
1.5 X 10" W < |Am|2| < 3.4 x 10" W 

0.36 < sin2^23 < 0.64 (2) 

at 90% CL, some important points remain unknown. We still lack information on the absolute 
neutrino mass scale uiq, only an upper bound mg ^ 0.2 — 0.7 eV [l^ exist. The neutrino 
mass pattern is not yet completely established: we do not know if nature prefers the normal 
(m3 > m2 > mi) or the inverted (m2 > nii > rris) mass hierarchy, where mi (777,3) is the 
mass of the neutrino state most (least) populated by the component. Moreover, we only 
have an upper limit on the mixing ang le ^13, sin2^i3 < 0.04, given by the CHOOZ reactor 



experiment 13|. 



The supernova density profile is such that a neutrino oscillation resonance in the 31-channel, 
involving Amg^ and 613, is bound to happen. Whether this happens for neutrinos or anti- 
neutrinos depends on the neutrino mass hierarchy (if normal or inverted). Besides, the reso- 
nance is very sensitive to the parameter ^13. The value of the hopping probability changes from 
zero to one when ^13 goes from 10"^ to 10"^. Clearly, supernova neutrinos provide an excellent 
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chance to determine two of the neutrino unknowns: the mass hierarchy and Oi^. This fact has 
been pointed out and explored by several authors 0, Q, 16]. 

However, the open questions on the dynamics of the SN explosion could, in principle, plague 
the determination of neutrino properties. One might wonder if it will be possible to disentangle 
the uncertainties of the supernova physics from the uncertainties on the neutrino parameters, 
and use experimental measurements of the SN neutrino fluxes to extract new information on 
both, the SN explosion mechanism and the neutrino oscillation parameters. 

In this article we analyze in detail the prospects for extracting the SN parameters as well 
as the neutrino oscillation parameters at three different types of next-generation detectors, 
from the measurements of neutrinos from the cooling phase of a galactic supernova. The most 
realistic next-generation experiments under present consideration are a megaton-scale water 
Cherenkov, a 100 kton liquid Argon and a 50 kton scintillation detector. We will study the 
performance of each of these detector types. In our analysis we will vary a total of seven 
parameters. Five are SN parameters: the average z/g, z/g and energies, respectively, (Ee), 
{Eg) and (E^); the ratio of the luminosities in x and e flavors, ^ (we assume the z/g and Ug 
luminosities to be equal) and finally the overall normalization of the fluxes, fixed by the total 
energy released (Efy) and the distance to the exploding star {D), Efy/D"^. The last two are 
neutrino oscillation parameters: the value of the angle ^13 and the neutrino mass hierarchy. 
We will simulate an observed set of data for given values of these seven parameters and use a 

method to henceforth construct confidence levels for the determination of these unknowns. 
In particular, we perform a comparison, highlighting strengths and weaknesses of each type 
of proposed experiment. In this paper we only consider the performance of the detectors 
regarding supernova neutrinos. However, when making a decision of which detectors should be 
build, their sensitivity to many other processes, including nucleon decay, solar and atmospheric 
neutrinos, beta-beams and super-beams as well as a neutrino factory, should naturally also be 
analyzed |17|. 

There has been a number of earlier works on supernova neutrinos, with many papers dis- 



cussing only the extraction of neutrino parameters. Refs. [18|, [lOj, however, investigate the 
possibility to get information on SN physics from the SN neutrinos using Super-Kamiokande 
and SNO detectors. These analyses take into account fewer supernova parameters and only 
some of the detection channels (considered as inseparable) we will consider here. Also, the 
simultaneous analysis of both normal and inverted hierarchy was not performed. The analysis 
in Ref. 20|] is very similar to the present study but concerns only a liquid Argon experiment. 
Bounds on neutrino masses, from the delayed time-of-flight as a function of the neutrino energy, 
obtainable with future large water Cherenkov as well as with a liquid scintillation experiment 
has also been discussed 211]. Other methods to extract information from a SN are: analyzing 
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23l . |24| . I25I . l26l | which may occur if the supernova neutrinos 



the earth matter effects 

traverse the Earth (mantle /core) before reaching the detector; studying the variation of par- 



n 



ticularly constructed variables [15|, |27|, |28|] , such as ratios of average e nergies, and recently the 
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33|. 



possibility of observing shock wave effects has attracted attention [15|, |29|, 

It should be noted that the determination of the SN parameters from the other studies, 
such as Earth matter effects are difficult, without prior knowledge of the value of 613. This 
is due to the fact that the Earth matter effects dependent on a combination of the hopping 
probability and the difference in unoscillated neutrino fluxes. Therefore, apparently a full 
analysis varying both neutrino and SN parameters is necessary for obtaining information on 
the SN parameters. In the present paper we will show to which maximal accuracy the SN 
parameters can be determined by the three types of experiments and how their determination 
depends on the values of the unknown neutrino parameters. Such a determination will likely 
be very helpful for the understanding of the physics of the core-collapse. 

In the next section we present the parameterized neutrino flux from the SN cooling phase. 
In Sec. mil we will discuss the analysis method and the three experimental setups we consider 
here. In Sec. |IV]we discuss our results and we devote the final section to our conclusions. 

II. THE NEUTRINO FLUX FROM THE COOLING PHASE 



In this analysis we will consider neutrinos emitted from the cooling phase of a type II 
supernova. This phase has the largest emission of energy with an expected total (time inte- 
grated) luminosity. Eh, of about 1 — 5 x 10^^ ergs. This luminosity is divided into all 6 flavors 
and we denote the individual contributions by Lj. We assume (as usual) that z/^, z/,-, P^, u^. are 
indistinguishable^ and denote them by the common index x. Thus, we have Eh = Le + Lg+A L^. 
We will furthermore assume that z/g and Ue luminosities are identical, i.e., Lg = Lg, which holds 
approximately in most SN simulation. We allow for a violation of luminosity equipartition by 
defining the parameter 

e = ^- (3) 

In general, simulations compute a value of ^ between 0.5 and 2 j^, 0]. 

We use the pinched Fermi-Dirac distribution for the energy-spectra of the neutrinos emitted 
from the supernova (the unoscillated flux) 



(4) 



^ We neglect the small difference in Ui and z/i, i = ^,t, fluxes originating from weak magnetism effects. 
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where is a normalization function. The average energies, (Ei), are hnearly related to the 
temperatures once the pinching parameter, rji, is fixed. As explained in the introduction the 
supernova physics strongly suggests that these average energies follow a hierarchy; (Ef.) < 
(Eg) < (Ex). The SN simulations favor < ?7e,?7g < 3 and rj^ < 2. In our simulations we take 
all ?7's to be zero, corresponding to pure Fermi-Dirac spectra. This is a conservative choice, 
since superposition of narrower energy spectra will be easier for the detectors to distinguish. 
In the case the average energies vary with time, the value of rj will be smaller, thus justifying 
our choice. Naturally, it would be even better if one could vary also the pinching parameters, 
but this is beyond our scope. 

The unoscillated Ui fiux at distance D from the supernova is given by 

(5) 

where all the luminosities are proportional to the total binding energy Ef,. It is worth remem- 
bering that one will normally not be able to see an optical counterpart to the SN (if the SN is 
inline with the galactic center it will be obscured by dust), leaving the distance D unknown. 
A crude estimate of the chance that an optical signal of the SN can be seen is only about one 
out of four. This should be compared to the fact that less than four supernova are expected 
in our galaxy per century. Therefore, in this article we will suppose that it is the combination 
Efj/D"^ that will be constrained by the SN neutrino detection. This can easily be translated to 
a constraint on the total emitted energy if the distance D can be independently determined. 
Thus, the unoscillated fiux of neutrinos from a supernova is parametrized by 5 variables: (Ee), 
{E,), {E,),^8.ndE,/D^ 

The fiux produced in the interior of the star will change its fiavor composition when traversing 
the outer parts of the star, due to neutrino fiavor oscillations in matter. As we assume that 
the fi and r fiuxes are identical, the oscillated fiux will only depend on the Ue and z/g survival 
probabilities. These probabilities strongly depend on the unknown neutrino parameters: 613 
and the mass hierarchy. We will now summarize these effects assuming that the matter density 
of the SN scales as p ~ r~^, which seems to be the most realistic density profile^. We will use 
the approximation 

p(r)=C -lO'^f^^Xg-cm-' , (6) 



r 

and we take the value C = 4. Due to the mass gap (Am2i/| Amigl ~ 1/30) and the smallness 
of ^13, the dynamics of the 3-z/ system can be factorized as two 2-z/ sub-systems: a high (H) 
one, driven by Amli and 613, and a low (L) one, driven by Amli and 612- Correspondently, 



^For others density profiles see 34 1. 
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FIG. 1: The z/g, Pee, and z/g, Pee, survival probabilities as a functions of sin^ ^13 for four different 
neutrino energies. We show the cases of inverted (left panel) and normal (right panel) hierarchy. We 
have used sin^ 6*12 = 0.3 and |Am|i| = 3.0 x lO'^ eV^. 

two resonances can happen as neutrinos travel through the SN. The H resonance will occur 
for neutrinos (antineutrinos) in the case of normal (inverted) hierarchy. In this case, there is a 
possible nonzero hopping probability Ph for crossing between effective mass eigenstates. The 
L resonance will occur for neutrinos^ and due to the values of the solar parameters given in ([1]) 
it will be always adiabatic (the level crossing probability ~ 0). Thus, we will not discuss 
this resonance further. 

The hopping probability Ph can be parametrized as 



H 



exp 



sin^ 9 



13 



1.08 ■ 10^ 



7\ 2/3 



\^rnl,\ 

10-3 



2/3 



^1/3 



(7) 



where E is in MeV and Arn^^ in units of eV^ ^ . It should be noticed that the value of C (which 
is uncertain by about a factor of 4) affects the translation from Ph to 613. Therefore, depending 
on how well the value of C can be determined, one might have to take this uncertainty into 
account when putting bounds on 61^. Knowing the type of the progenitor star would be a help 
in this case. 

The z/g and Ue survival probabilities. Pee and Pge, respectively, are approximated for the 



■l-Since SNO's demonstration that the solar CC/NC ratio is less than | we know that Arn^i > 0. 

''^ Strictly speaking the Am^ in Eq. ^ should be Am^2 for the normal hierarchy and Ani^i for the inverted 



one. However since jAml^l 



jAmljl, this is not relevant here. 
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normal hierarchy by 

Pee ^ PH\Ue2\'+il-PH)\Ue3\', (8) 
Pee ^ \Uel\^ , (9) 

where Ua,i ,a = e,fi,T ,i = 1,2,3 are the Maki-Nakagawa-Sakata neutrino mixing matrix 
elements and we have used here the standard parameterization. Using ([2]) we find |f/eiP — 0.7, 
|t/e2p — 0.3 and |f/e3p < 10~^. And for the inverted hierarchy 

Pee ^ \Ue2\^ (10) 
P,, ^ PnlUell' + (1 - PH)\Ue3\' • (11) 

These survival probabilities as a function of sin^ ^13 are shown in Figure [H for various neutrino 
energies and the two types of hierarchies. We will neglect collective fiavor transformation of the 
neutrinos, which in a detailed analysis should be added to the conventional MSW conversion 



in the supernova environment 35j 



The final fiuxes arriving at Earth is simply given by 

F.^ = FOp,e + i^°(l-Pee), (12) 

Fp, = F°Pg, + F°(l-P,,), (13) 

F,^ + F,^ = F° (1 - Pee) + F° (1 + Pee), (14) 
Pp, + Pp. = Flil- P-e,) +Fl{l + P-e^) ■ (15) 

Notice that in the case of luminosity equipartition and degeneracy in average energies, the fiux is 
independent of the survival probabilities and thus also independent of the neutrino parameters. 

When calculating the oscillation probabilities we have neglected the Earth matter effects [l^ . 
These effects can be precisely calculated and thus will not affect our conclusions much, if the 
direction of the supernova is known. Even in the case that an optical counterpart of the SN 
can not be seen, the direction of the supernova can be determined rather well from the neutrino 
fiux, if the detector (or another existing detector than the one being analyzed) can measure the 
elastic scattering of neutrinos on electrons. This detection channel is highly forward peaked and 
in the case of Hyper-Kamiokande the direction can be inferred to within ~ 1° [3^. Therefore, 
we will not take into account the Earth matter effects in our calculations. However, in section 
IIVI we will briefly discuss the possibility of extracting further information about the hierarchy 
by detecting Earth matter effects. 
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III. THE ANALYSIS METHOD AND THE DETECTORS 



In this section we will describe our method for studying the sensitivity of the detectors to 
the parameters under investigation of neutrinos from a nearby supernova. While we are rather 
optimistic in our choice of detector performances we will be taking somewhat difficult choices 
for the input SN parameters. 

As explained earlier, the neutrino fluxes arriving at the detectors will depend on 7 parameters 
(5 SN dynamics parameters, and 2 neutrino physics parameters). For each experiment, we 
simulate the expected number of neutrino events at each observable mode, for a fixed set of 
these parameters. The artificially generated data will be our input data (the imagined true 
values). We then construct a = X^((-^e), {Eg), {Ex) , ^, Eb/ D"^ , siii^ 613, sign{Aml-^^)) function 
in order to fit these unknown parameters to the input data in the usual way. We only consider 
statistical uncertainties in our x^- This allows to compare the maximal attainable sensitivity 
for each detector type. We compute the allowed regions for each pair of unknown parameters 
to estimate the experimental sensitivity to them by marginalizing with respect to the other 
4 parameters, for a fixed hierarchy. Since we construct the confidence level region for each 
hierarchy separately, our graphs will still be useful, in the case that the neutrino mass hierarchy 
is determined before a SN observation. Moreover, this allows us to conclude whether or not 
the hierarchy can be established. 

The artificially generated data are constructed simply by calculating the theoretical expec- 
tation at the chosen set of input parameters. We have tested that this gives similar results as 
when a more realistic data set is used. We have constructed a Gaussian distributed data set, 
by choosing randomly a point from a Gaussian distribution centered around the expectation 
for the bin and of width equal to its square-root. There is no difference in the results from 
the Gaussian distributed data set and the data set given by the theoretical expectation, be- 
sides that the Gaussian data set might give a fluctuation of the allowed regions away from the 
central value. Generating the observed data set by the theoretical expectation allows for easy 
independent reconstruction and comparison and we therefore prefer to use this simpler method 
for generating the data. 

We will in this analysis use the following parameters space when varying the SN parameters: 

(Ee) e [9 - 15] MeV , (Eg) G [12 - 17] MeV , (E,) e [15 - 30] MeV , 

^ e [0.5 - 2.0] , Eb/D^ e (2.0 - 4.0) X 10^^ ergs/kpc^ . (16) 

This parameter space is roughly what is expected from SN simulations . In any case if an 
analysis of SN neutrino observation cannot determine the parameters within the area suggested 
in (|T6i) . the SN simulations will probably do a better job. Naturally, the confidence levels and 
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the accuracy with which each parameter can be measured, will depend on the input data. For 
instance, in the case when (Ef.) = (Eg) = (E^) and ^ = 1 the neutrino flux becomes independent 
of the survival probabilities and thus of the neutrino oscillation parameters. Therefore, in this 
special limit (luckily not favored from SN simulations) the neutrino parameters cannot be 
deduced from the observation of SN neutrinos. We will discuss how the neutrino parameter 
determination depend on the strength of the hierarchy of the average energies. 





(^e) (MeV) 


{Ee) (MeV) 


{E,) (MeV) 




Eb/D"^ (ergs/kpc^) 


point 1 


12 


15 


18.0 


1.50 


3.0 xlO^i 


point 2 


12 


15 


18.0 


0.75 


3.0 xlO^^ 


point 3 


12 


15 


16.5 


1.50 


3.0 xlO^^ 



TABLE I: Definition of the input SN parameters at some reference points. 



In Table |T] we define three points in the SN parameter space, that we will use as reference 
points. These points are situated within the expectation of the SN simulation given in equation 
[TBI Most parameters are chosen at their central value. However, the value of ^ has been chosen 
at two point slightly away from its central value (being 1), so as to see its impact on the 
accuracy with which the SN and neutrino parameters can be determined. Also the values of 
(E^) is chosen at the lower range of its expectations. This has been done in order to see until 
which lower value a meaningful determination of the neutrino parameters can be expected. It is 
important to notice that increasing the value of (E^) will increase the sensitivity of the detectors 
to the neutrino parameters as well as to the supernova parameters. First of all a higher value 
of (Ex) will give a larger number of events, since the detection cross-sections increases with 
energy. Also the increased hierarchy between (E^) and (Eg) will make the impact of the value 
of the survival probability larger. For points 1 and 2 we study a rather conservative case of a 
hierarchy of only 19% difference between (Eg) and (E^). For point 3 the hierarchy is of only 
10%, in which case it will become harder to determine the neutrino parameters. Although, such 
a weak hierarchy between the average energies seems rather unlikely it cannot yet be excluded. 

The solar neutrino mixing angle and the atmospheric mass-square difference are fixed to 
values in the allowed region given in ([2]): 

= 0.575 rad , |Am^J = 3.0 x 10"^ eV^ . (17) 

The atmospheric mixing angle does not enter into the calculation as and Ut enter on the same 
footing both in the SN and in the detector (it introduces a unobservable rotation). Furthermore, 
as we do not take into account Earth matter effects, we have no dependence on the solar mass- 
square difference. When calculating the confidence levels we will just study the extreme cases 
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for the true value of 613. Meaning that we will only investigate two values of Ph', zero and one. 
In this case we have three scenarios for the input neutrino parameters: 



1. Scenario iO: Inverted hierarchy and Ph — (large ^13); 

2. Scenario nO: Normal hierarchy and Ph — (large 6*13); 

3. Scenario al: Any hierarchy and Ph — 1 (small ^13). 

In the case of Ph — 1 the inverted and normal hierarchy are identical. In the above scenarios 
large ^13 means a values corresponding to sin^ 9i-^ = 10~^ and small means a value corresponding 
to sin^6'i3 = lO"*^ (see figure [I]). 

We will analyze three different types of next-generation experiments, namely: 

• Water Cherenkov (WaterC); 

• Liquid Argon (LAr); 

• Scintillation. 

Each of these will be discussed in detail in the following subsections. The WaterC and Scintilla- 
tion experiments are much more sensitive to anti-neutrinos than neutrinos, due to the dominant 
inverse beta decay detection. The dominant detection channel for a LAr detector on the other 
hand is charged current z/g interactions on Argon, making it more sensitive to neutrinos. 

Before discussing the details of each detector we will note some common features. First of all, 
obviously having sensitivity to more than one combination of neutrino fluxes will be essential 
to pin down the SN and neutrino parameters. For each experiment we will make contours for 
two scenarios: a pessimistic one, assuming that only the dominant channel can be used, and 
an optimistic one, where we assume several channels can be separated by the detector. This 
will illustrate the necessity of having sensitivities to several channels with different sensitivity 
to z/g/ Uf. and fluxes. In particular, a NC channel along with a CC channel will complement 
each other well. A true NC channel is independent of the neutrino parameters, as is easily seen, 
as we have 



which is given purely in terms of the original fluxes. However, it should be remembered that 
eg. elastic scattering on electrons is a combination of NC and CC interactions for the electron 



iVNC OC / {F^,0-NC + Pp,(^NC + IF^^aNC + 2Fp^(TArc) dE , 



(18) 




(19) 
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neutrino and anti-neutrino species. Thus for this channel cxp^ 7^ a^^ 7^ a^^ and the number of 
events have a shght variation with the oscillation parameters. 

We neglect any time-dependence of the neutrino fluxes and just look at the energy spectrum, 
or in the case that the energy cannot be measured, the total number of events integrated 
over time. The time-dependence of the energy spectra can be monitored by the experiments 
we study. Therefore, it will be possible to get a feeling on how well the time-independence 
assumption works. The SN simulations suggest that a mild steady increase in the average 
energies as a function of time will occur. This suggest that indeed a pinched Fermi-Dirac 
spectrum, can still be used, but with a broader spectrum and thus a smaller value of rji. Also 
a steady decrease in the luminosities is expected, not influencing our method. 

We will also assume that all detection efficiencies are 100% above the threshold. This is 
a very optimistic assumption, but the efficiencies for these future experiments are presently 
unknown, and at least in this manner we treat all three experiments on the same footing. 
Also, we will neglect the energy resolution, but this is partly compensated by the use of wide 
bins in energy. We will use energy bins of 10 MeV, unless otherwise specified. The exact 
energy resolution, which is also unknown, can become important in some cases where there are 
degeneracies between certain parameters. However, for most of the parameter space the exact 
energy resolution is not very important and will not change the general result". 



A. Analysis of a Water Cherenkov Detector 



The water Cherenkov detectors have proven very successful, with the Super-K collaboration 
being first at announcing extremel y c ompelling evidence of atmospheric neutrino oscillation 
in 1998. The Hyper-Kamiokande [371] (Hyper-K) detector is being proposed to replace the 
current Super-Kamiokande experiment and will have a total mass of about one megaton. Other 
proposed and more or less identical detectors, with the only difference being their location 
and the exact mass, are the (American) UNO detector 38| and the (European) MEMPHYS 
detector 39|]. All these detectors will of course have similar sensitivities. We assume the fiducial 



volume of the WaterC to be 540 kton, which is the expectation for the Hyper-Kamiokande 37] 



detector. Earlier works on the subject can be found in |18l . |19| although these references only 
take into account one detection channel. As we will prove, the possibility to measure neutral 
current and charged current on oxygen as well as elastic scattering on electrons, will greatly 
improve the sensitivity of a WaterC detector. 



" The energy resolution is, however, important for the determination of Earth matter effects for supernova 
neutrinos from a single detector 
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FIG. 2: The cross-section for the various channels for neutrino detection in a WaterC detector. Here 
ES stands for elastic scattering on electrons. 



We will take into account four different channels. Two channels are CC reactions that will 
provide spectral information for Ve and Ve fluxes. One channel is a NC reaction and thus 
sensitive to all neutrino flavors, providing information on the total neutrino flux. The ELAS 
channel is, as discussed earlier, also mainly an NC reaction, but with a small CC contributions. 
Below we list the four channels. 

1. The inverse beta decay (IB) for detection of z/g 

z>e + P ^ e+ + n . (20) 

d to be 5 



This is by far the dominant detection channel and we will assume the thresho^ 
MeV. The cross-section is well-known and we use the calculation given in Ref. 



40| 



2. The absorption of Ve and Ve on oxygen by CC interactions (CC-0) 

+^^^0 ^ e± + X . (21) 



4l| and the threshold for detection is about 15 



The cross-sections are taken from Ref. 
MeV. 

3. The elastic scattering (ELAS) on electrons 

Vi + e^Ui + e , (22) 
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is possible for all types of neutrinos, although the cross-sections for u^. and z/g are slightly 
higher due to the additional CC contribution. The ELAS channel is easily separated as 
these events are strongly forward peaked. We have set the detection threshold at 7 MeV. 

An interesting channel for observation is the excitation of oxygen by NC interactions, 
followed by a decay chain with emission of a detectable mono-energetic photon, as first 



discussed in Ref. j42|] 



Vi+^^O^Vi + -i + X . (23) 

Hyper-K can detect photons with an energy greater than about 5 MeV. The excited ^^O 
atom decay to either ^^O* or ^^N*, which then emits a photon with the energy in the range 
5-10 MeV. As the photon has a well defined energy, these events can be easily separated 
from the other detection channels. We will refer to this channel as NC-0, even though 
it should be remember that it only includes those partial decay chains which give rise to 
a detectable mono-energetic photon. This channel differs from the others as the energy 
spectrum cannot be measured. Henceforth we only use the total number of events. 

We will for the main parts of this article assume that the four channels can be separated, and 
make remarks on the case where only the inverse beta decay channel is detectable. The NC-0 



and ELAS channels should be easily separated from the others. As discussed in Ref. 3l| , the IB 
events should be slightly forward peaked, whereas the CC-0 events should be slightly backward 
peaked and this gives an opportunity to distinguish these events. Moreover, by addition of small 



amounts of gadolinium 43|] in the WaterC detector, the capture of neutrons is possible and this 
would assure the separation of the IB and CC-0 channels. This is also the reason that we 
have put the detection threshold for the IB channel as low as 5 MeV. Presumably even a lower 
threshold can be achieved with a gadolinium enriched WaterC detector. For the IB, CC-0 and 
ELAS channels we calculate the energy spectrum and we use 10 bins with a width of 10 MeV 
and the first bin starting at 5 MeV. In Figure [2] we show the cross-sections involved. 

In Table im we show the total number of events we calculate for the SN parameters at point 
1. The dominant channel for Hyper-K is the inverse beta decay channel, which is only sensitive 
to the Ve flux arriving at the Earth. However, the NC-0, CC-0 and ELAS channel is very 
important for having sensitivity to other combination of fluxes. Moreover, the ELAS channel 
gives the opportunity to determine the direction to the supernova and thus an early warning to 
astronomers will be possible. From Table [TTl it is expected that if the true neutrino parameters 
are consistent with the inverted hierarchy and large ^13 then Hyper-K will be able to determine 
this scenario with a very high confidence level. But, even the sub-dominant channel allows for 
a determination of the normal hierarchy if the angle 6*13 is sufficiently large. In section IIVI we 
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Expected Number of Events in a 540 kton water Cherenkov detector 

{Ee) = 12 MeV, {Ee) = 15 MeV, {E^) = 18 MeV, = 1.5 Lg 


Reaction 


Any hierarchy 
Small 


Inverted hierarchy 
Large 6*13 


Normal hierarchy 
Large (9i3 


Inverse beta decay 


1.4 X 10^ 


2.1 X 10^ 


1.4 X 10^ 


CC on oxygen 


7.7 X 10^ 


10.7 X 10^ 


9.5 X 10^ 


v^ + e^Vx + e 


8.4 X 10^ 


8.7 X 10^ 


8.8 X 10^ 


NC on oxygen 


3.5 X 10^ 


3.5 X 10^ 


3.5 X 10^ 


TOTAL 


1.63 X 10^ 


2.34 X 10^ 


1.65 X 10^ 



TABLE IL Expected number of neutrino events in a 540 kton water Cherenkov detector for the SN 
parameters at point 1 and each of the cases al, iO and nO for the neutrino parameters. 

will comment on how this dependent on the SN parameters (luminosities and average energies), 
by using our reference points. 



B. Analysis of a LAr Experiment 

Next, we will look at the possibility to determine the supernova and neutrino parameters 
at a future Liquid Argon experiment. The are various LAr experiments proposed. The Icarus 
detector at CNGS is expected to have a 3 kton final version^ and already a 300 ton detector 
is running. Moreover, the LANNDD 4J], the GLACIER 45| and the Flare 46| detectors are 
being discussed as possible 
can be found in 



20 



28 
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'uture detectors in the 100 kton size. Earlier works on the subject 



48 



We will take into account the following channels: 
1. Detection of Ve through CC- interaction: 

Ve Ai^e- + A' + nN (Ethr = 1.5MeV) 



(24) 



where nN represent emitted nucleons or other debris (like a-particles etc.) and A' is the 
leftover nucleus. The CC cross-section on Argon are taken from Ref. 481. l49|. 



2. Detection of v^. through CC-interaction: 

z7, Ar 
having a threshold about 7.5 MeV. 



A' + nN , 



(25) 



t Unfortunately, ICARUS might be interrupted due to cancellation of its funding. 
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1 00 



E (MeV) 

FIG. 3: The cross-section for the various channels of neutrino detection in LAr. 



3. The elastic scattering (ELAS) on electrons 



iy^ + e. 



(26) 



is possible for all flavors of neutrinos. Again we have taken the threshold to be 7 MeV. 
4. The scattering on Argon of any type of neutrino through NC-interaction: 



Ui At ^Ui + A' + nN. 



(27) 



This channel has no sensitivity to the energy. The NC cross-section on Argon are taken 
from Ref. 
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49|. 



As in Ref. 



20| it is assumed that one can separate all four channels. However, we will 



also make contours for the 'worse case' scenario were only detection in the z/g CC channel is 
available. For the z/g and z/g CC reactions and the z/j NC reactions, the energy and time-delay of 
the photons emitted from the de-excitation of respectively K, CI and Ar can be used to classify 
the type of event. For the ELAS events no photons will be present. The cross-sections are 
shown in Figure [3l and it should be noted that we do not take into account in our calculations 
their uncertainties. At the moment there are no experimental confirmation of the theoretically 
calculated cross-sections. But hopefully, in the case that a large scale LAr detector will be 
realized, the cross-sections will already have been experimentally measured (eg. by ICARUS). 
For the channels with measurable energy spectra we again use 10 energy bins of 10 MeV each, 
the first bin starting at 5 MeV. 
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Expected Number of Events in a 100 kton LAr detector 

{Ee) = 12 MeV, (Eg) = 15 MeV, {E^) = 18 MeV, L.^ = 1.5Le 


Reaction 


Any hierarchy 
Smah 


Inverted hierarchy 
Large 6*13 


Normal hierarchy 
Large (9i3 


fe CC on Argon 


1.4 X 10^ 


1.4 X 10"^ 


1.7 X 10^ 


i/g CC on Argon 


4.2 X 10^ 


7.9 X 10^ 


4.2 X 10^ 


ELAS 


1.2 X 10^ 


1.3 X 10^ 


1.3 X 10^ 


NC on Argon 


1.3 X 10^ 


1.3 X 10"^ 


1.3 X 10^ 


TOTAL 


2.80 X 10^ 


2.84 X 10"^ 


3.16 X 10^ 



TABLE IIL Expected number of neutrino events at a 100 kton LAr detector for the SN parameters 
at point 1 and each of the cases al, iO and nO for the neutrino parameters. 



In Table IIIII we show the total number of events for the SN parameters at point 1 for each 
detection channel. The dominant channel is detection of z/g by the charged current interaction 
on Argon. Also the NC and ELAS channels have a fairly large number of events, whereas the 
sensitivity to the flux is rather weak. 



C. Analysis of a Scintillation Detector 



Finally, we will examine the proposal for the 50 kton Low Energy Neutrino Astronomy 
(LENA) 50|, l5l| liquid scintillation detector. We assume that LENA will be filled with pure 
PXE (CigHig). If another oil will be used the carbon to proton ratio may change and thus 
the results will change slightly. A discussion on SN neutrinos and scintillator detectors can be 
found in 5ll, l52] , although without an explicit calculation of the accuracy of the determination 



of the parameters. 

We will exploit six v detection channels in LENA (three are CC reactions, two are NC 
reactions and the last is the ELAS channel), these are listed below: 



1. The inverse beta decay for detection of 



P 



n 



(2J 



The threshold energy for this reaction is 1.8 MeV. Again we take the IB cross-section 
from |40|]. 

2. The CC capture of on ^^C 



V, V C ^ ^^B + 



(29) 



12b c + e- + z^e 
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10' 



o 



to lo-^P 



u.+C (CC) 




1 00 



E (MeV) 

FIG. 4: The cross-section for the various channels of neutrino detection in LENA. 

The threshold Ue energy for the capture on ^^C is 14.39 MeV. 

3. The CC capture of on ^^C 

c + e+ + i/e 

The threshold neutrino energy for capture on ^^C is 17.34 MeV. 

4. Elastic scattering on protons 

Ui+p ^ Ui+p (32) 

This process might in some areas of parameters space even give a larger number of 
events than the IB process, due to the factor of six, originating from the number of 
neutrino and anti-neutrino species. The cross-sections is taken from 



(30) 
(31) 



where we have 



used the approximation of equal cross-section for neutrinos and anti-neutrinos. We have 
implemented a neutrino threshold energy of 25 MeV, corresponding roughly to a cut of 
0.2 MeV in electron equivalent energy. 



5. NC scattering on C: 



, 12 



(33) 



I2n* 



,12 



C + 7 (15.11 MeV) 



18 



Expected Number of Events in a 50 kton Scintillation Detector 

{Ee) = 12 MeV, (^e) = 15 MeV, {E.^) = 18 MeV, L.^ = 1.5Le 


Reaction 


Any hierarchy 
Small 


Inverted hierarchy 
Large 


Normal hierarchy 
Large (9i3 


Inverse beta decay 


1.0 X 10^ 


1.5 X 10"^ 


1.0 X 10^ 


z/g CC on carbon 


6.0 X 10^ 


1.1 X 10^ 


6.0 X 10^ 


Ue CC on carbon 


1.0 X lO-'^ 


1.0 X 10^ 


1.4 X 10^ 


Vi NC on proton 


9.9 X 10^ 


9.9 X 10^ 


9.9 X 10^ 


Vi + e ^ Vi + e 


7.9 X 10^ 


8.2 X 10^ 


8.2 X 10^ 


Ui NC on carbon 


1.4 X 10^ 


1.4 X 10^ 


1.4 X 10^ 


TOTAL 


2.39 X 10^ 


2.93 X 10"^ 


2.43 X 10"^ 



TABLE IV: Expected number of neutrino events in a 50 kton scintillation detector for the SN param- 
eters at point 1 and each of the cases al, iO and nO for the neutrino parameters. 



The emission of a mono-energetic photon, makes this channel easily separated from the 



others. This cross-section is taken from S^]- Since the emitted photon carries no infor- 



mation about the neutrino energy, this is the only channel for which LENA will have 
no energy information. Therefore, we will only use the total number of events from this 
process. 

6. EL AS on electrons 

z/j + e — > z/j + e , (34) 
which has been discussed in earlier sections. 

Here again we take two approaches. The conservative one, where we consider only the 
inverse beta decay channel and the optimistic one, where we assume that all channels can be 
distinguished from each other. In principle, one can hope that, due to the distinctive signatures 
of the above discussed channels they can be separated. The most doubtful discrimination is 
between the z/g and Ve CC reactions on carbon. It might be possible to separate these by using 
the delayed coincidence of the (3^ / [3~ decays with the primarily produced electron/positron 
and the knowledge of the average lifetimes of the produced unstable nuclei. For the channels 
with sensitivity to the neutrino energy spectra (all except the NC scattering on carbon) we 
again use 10 energy bins of 10 MeV each, the first bin starting at 5 MeV. Moreover, due to 
the very fine energy resolution expected for LENA, with a threshold of order 200 keV, we also 
include an extra low-energy bin, with the events originating from neutrino energies below 5 
MeV, for the IB and ELAS channels. 

The dominant channels for LENA are the inverse beta decay channel (for z/g) and the z/j NC 
scattering on protons. In Table [IV]we show the total number of events for the SN parameters at 
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point 1 for each observable neutrino channel in LENA. The detector will have a large sensitivity 
to both z/g and the total neutrino flux and even a reasonable sensitivity to the v^. flux. 

IV. RESULTS 

In this section we will study two scenarios for the detection capability for each detector: the 
case where only the main CC channel is available and the case where all channels discussed in 
section IIIII are available and moreover separable. As detecting only the main CC channel is a 
rather pessimistic scenario, we include in appendix A a discussion of the detector performance 
in a few other cases with a more realistic, but non-optimal, detector-setup. The main reasons 



for studying the main CC channel only case is to compare to other studies (Refs. 18|, ll9|]) and 
also to illustrate the importance of having both NC and CC channels. Furthermore, this is the 
absolute worst case scenario and with the two contours considered in this section we give the 
span of possible allowed regions. 

In Figs. [5], [6] and [7] we show the Sa CL allowed regions obtained for a WaterC detector in 
the three limiting cases al, iO and nO, discussed in Sec. IIIII The corresponding regions for the 
LAr and the Scintillation detectors are shown in respectively Figs. [HI E] and [ID] and Figs. [HI [T2] 
and [131 

In the case al the hierarchy cannot be determined, as both hierarchies produce the same 
neutrino fluxes. All three experiments present a good sensitivity to ^13 as can be seen from 
Figs. [SI [51 and [TT] and will provide an upper limit on sin^ 6*13 of about 1-2 xlO~^, if there are 
several channels available. In the case that only the main channel is available (IB or z^eCC) a 
degeneracy in 6*13 and the hierarchy occurs and 6'i3 cannot be determined (unless the hierarchy 
has already been established by another experiment). This independence of ^13 for the case of 
normal hierarchy and the WaterC and Scintillator detectors, is evident from FigUl as the IB 
channel is only sensitive to Pgg, which is almost constant as a function of 6*13. Similarly, the 
LAr detector is not sensitive to ^13 in the case that only the i^eCC channel is available and the 
hierarchy is inverted. The anti-neutrino flux arriving at Earth for the case al (and nO) consist 
of roughly 70% with a temperature of 15 MeV and 30% with a temperature of 18 MeV (for 
point 1). For the Scintillator detector the statistics for the only IB case is such that a pure flux 
with temperature of about 16.5 MeV can also flt this data. Therefore, the inverted hierarchy 
is allowed even for large values of ^13. However, for a small interval around sin^ ^13 ~ 10~^ and 
low values of (-Eg) it is not possible to flt the energy spectrum, giving rise to islands in the 
contours as seen on FigHH (the exactly same contours are found in Fig JT3l) . If sin^ 6*13 turns out 
to be in the region where Pu — 1, then there is no possibility that it can be determined by any 

In this case SN information 



of the laboratory experiments currently proposed 55 



56 



57 
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FIG. 5: Sensitivity of a 540 kton WaterC detector assuming true SN parameters as in point 1 of 
Table m for any hierarchy and sin^(0i3) = 10"^ {Ph — !)• We show 3cr CL contours (2 dof) using: 
all 4 channels (IB + CC-0 + ELAS+ NC-0) and normal hierarchy (NH) marked by the dark (blue) 
horizontally hatched area; only the IB channel and normal hierarchy marked by the dark (blue) dashed 
line; all 4 channels and inverted hierarchy (IH) marked by the light (cyan) diagonally hatched area; 
only the IB channel and inverted hierarchy marked by the light (cyan) dashed line. 
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FIG. 6: Same as Fig. [5] but for the inverted hierarchy and sin2(0i3) = 10-3 {Ph ~ 0). The normal 
hierarchy is ruled out by more than 5(T, having a global Xmin — case of all four channels 

being present. 
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FIG. 7: Same as Fig. [5] but for the normal hierarchy and sin^(0i3) = 10~^ {Ph — 0). The inverted 
hierarchy is ruled out by more than 5(T, having a global Xmin ~ case of all four channels 

being present. 
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FIG. 8: Sensitivity of a 100 kton LAr detector assuming true SN parameters as in point 1 of Table 
m for any hierarchy and sin^(0i3) = 10~^ {Ph — !)• We show 3(7 CL contours (2 dof) using: all 4 
channels (z^eCC+ z^gCC + ELAS + NC) and normal hierarchy marked by the dark (blue) horizontally 
hatched area; only the z^eCC channel and normal hierarchy marked by the dark (blue) dashed line; all 
4 channels and inverted hierarchy marked by the light (cyan) diagonally hatched area; only the fgCC 
channel and inverted hierarchy marked by the light (cyan) dashed line. 
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FIG. 9: Same as Fig. [8] but for the inverted hierarchy and sin2(0i3) = 10-3 {Ph ~ 0). The normal 
hierarchy is ruled out by more than Aa, having a global Xmm ~ in the case of all four channels 
being present. 
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FIG. 10: Same as Fig. [8] but for the normal hierarchy and sm^iOvi) = 10-3 {Ph ~ 0). The inverted 
hierarchy is ruled out by more than 5(T, having a global Xmin ~ case of all four channels 

being present. 



26 









— ^ 


* 


^^-\\ '\ 





10 



10 



10 



10 



sln'(0,3) 




sln'(0,3) 




50 kt Scintillation (a1) 

1^ Input (any hierarchy assumed) 

— IH 3(7 

IH 3 CT ONLY IB 

— NH3a 

■■■■ NH3CT0NLYIB 



E,/D' (10'' Ergs/kpc') 



FIG. 11: Sensitivity of a 50 kton Scintillation detector assuming true SN parameters as in point 1 of 
Table m for any hierarchy and sin^(0i3) = 10"^ {Ph — 1)- We show 3a CL contours (2 dof) using: all 
6 channels (IB + feCC+ z/gCC + z/ — p + NC + ELAS) and normal hierarchy marked by the dark 
(blue) horizontally hatched area; only the IB channel and normal hierarchy marked by the dark (blue) 
dashed line; all 4 channels and inverted hierarchy marked by the light diagonally (cyan) hatched area; 
only the IB channel and inverted hierarchy marked by the light (cyan) dashed line. 
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FIG. 12: Same as Fig. [11] but for the inverted hierarchy and sin^(0i3) = 10~^ {Ph — 0). The normal 
hierarchy is ruled out by more than 5cj, having a global Xmin ~ ^'^^ case of all six channels 

being present. 
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FIG. 13: Same as Fig. [TT]but for the normal hierarchy and sin2(0i3) = 10-3 (p^ ~ 0). The inverted 
hierarchy is ruled out by more than 5a, having a global Xmin ~ case of all six channels being 

present. 



29 



will be extremely valuable. If several channels can be included in the analysis, the value of the 
SN parameters {E^)-, (Eg) (except for the LAr detector) and Eh/D"^ can be determined to a few 
% level in case al. On the other hand, ^ and {E^,) will be much less constrained by data (for 
point 1). 

If sin^ is large, i. e. in the region where Ph ~ (cases iO and nO), all three experiments can 
determine the neutrino mass hierarchy. From Figure [6] and Figure [7] we clearly see that if ^13 is 
large, the neutrino mass hierarchy can be determined with a high confidence level independent 
of the true hierarchy, in a WaterC detector. Also the LAr and the Scintillation detectors can 
establish the hierarchy in the large 6*13 region, as shown in Figs. [9lfT0l and [T2lfT3| although with 
less significance. It should be noted that, as we will discuss below, the establishment of the 
hierarchy depends on the parameter space assumed for the SN parameters, and in particular 
whether (Eg) is allow to undertake the same value as (Ex). 

In the case iO, the WaterC detector can determine the hierarchy and will give a lower limit 
on sin^6'i3, even if only the IB channel is available. For the LAr detector and the scenario 
with only the z/gCC channel being present, we see that the normal hierarchy can explain the 
iO scenario. This is again evident from Fig{T] as in this case the information contained in the 
measurement is that Pgg — 0.3, which occur for the normal hierarchy and small values of ^13 and 
for inverted hierarchy and any value of 613. For the only IB case and the Scintillator detector, 
the measurement contain the information of the Ue flux with a temperature corresponding to 
the value of (E^). This can also be fitted with a normal hierarchy and raising the value of (E^) 
(giving about 30% of the flux) and (Eg) (giving about 70% of the flux) such as to simulate 
a pure z/g flux of 18MeV. The statistics of the WaterC detector is good enough to reject this 
case for the SN parameter space given in EqIT6l If all channels are considered in the data 
analysis, the LAr and Scintillation detectors can produce a similar constraint on sin^ ^13 as the 
WaterC detector and the lower limit is given by 1-2 xlO~^. For scenario iO the z/g flux arriving 
at the Earth is identical to the original flux and thus in this case one looses much of the 
sensitivity to (Eg). But, a slightly better sensitivity to (E^) is obtained as compared to the al 
case, where a part of Fp^ is the original z/g flux. The LAr detector can measure (E^) (even if 
only U(.CC events are used) and E^/D'^ to about 5 %, but will not be very sensitive to the other 
SN parameters. Also, the LENA-type detector, can determine (E^) and Eh/D"^ to a few % 
with the help of all channels. The WaterC detector will be very sensitive to all SN parameters, 
except for (Eg), if all channels contribute. 

In the case nO, even if only the main channels are available, (E^) can be determined to a 
few % by the WaterC detector (see FigJT]). In this scenario there is basically no sensitivity to 
{Eg), as the flux is identical to the original flux. But, all the others SN parameters can 
be accessed with very good precision by the WaterC detector, if all channels take part in the 
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analysis. The LAr setup can measure {E^) to a few % and {E^) and Eh/D"^ to about 10 % 
(see FigJTOl). A LENA- type detector will have a sensitivity of about 5 % for (E^), (Eg) and 
Eiy/D"^ and about of 20% for ^. Furthermore, if all channels are considered in the data analysis, 
the LAr and Scintillation detectors can again provide similar bounds on sin^ 6'i3 as the WaterC 
detector, giving a lower limit of 0.7-2 xlO~^. If only the IB channel is available not even the 
WaterC detector can determine the hierarchy. As seen from Figil] the value of Pgg of roughly 
0.7 is also found for the inverted hierarchy and small values of Oi^. 

Let us shortly compare the sensitivity to the ^13 angle to that of other proposed future 
experiments. One expects the reactor experiment Double Chooz to reach down to sin^ ^13 ~ 
5 X 10^'^, and the next generation ^13 reactor experiment Daya-Bay to reach sin^ ^13 ~ 2.5 x 
10~^ 551]. The proposed novel technique which exploits the recoilless resonant absorption of Ue 
to measure sin^ ^13 in a short baseline experiment, may be able to reach similar sensitivity 56|. 
So to reach sensitivities to sin^ 9is ^ 10""^ before a new galactic SN observation one would 



probably need beta-beams 57| or neutrino factories 58[]. Correspondingly, the determination 
of a lower bound on of order 10"^ from a SN observation for the cases iO and nO, can be of 
great importance. 

In the following we shortly compare to the work in Ref. 19j, where only the inverse beta 
decay channel at a WaterC detector was analyzed. In this paper only the case al is studied 
and the 6*13 angle is not varied. This seems unlikely as there are no known experiment that 
can restrict 6*13 to be smaller than 10~®, which would be necessary for getting independence of 
the exact value of the CHOOZ angle. Therefore, it is erroneously concluded in Ref. 19| that 
the value of ^ can be well determined at HyperK with only the IB channel available. As can 
be seen from Figj5] there is a degeneracy between ^13 and C,, only broken by the addition of 
NC channels. This indetermination of ^ in turn influence the measurement of E^, as there is a 
degeneracy between Ef, and ^ (see FigjS] panel 9). Therefore, overall the allowed regions found 
in Ref. [3] are too restrictive as compared to our contours for the case of only IB. 

Next we compare to the work in Ref. 18|], where the WaterC detector with the inverse beta 
channel along with the CC-0 was studied, although these were assumed inseparable. In addition 
the SNO detector was also analyzed, but as the fiducial mass is much smaller, this experiment 
does not add significantly to the statistics. The main sensitivity in Ref. 18[ is to the z/g flux 
with only a minor sensitivity to the Ug flux. In this paper it is concluded that the HyperK 
detector cannot determine the normal hierarchy even for large 613. Clearly, we don't agree 
with this statement, as with the addition of the sub-dominant channels this is indeed possible. 
However, we see that if only the IB channel is available this conclusion is true. But, it must be 
remembered that the only-IB scenario is not a realistic one. We also note that for small the 
very precise determination of (Eg) claimed in Ref. 18[ is only possible for the case that normal 
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point 2 






point 3 






iO 


nO 


al 


iO 


nO 


al 


{Ee) IH 


11.0-13.2 


— 


11.1-13.4 


10.0-14.1 


— 


< 14.2 


(Ee) NH 


— 


10.8-13.9 


10.5-13.0 


< 13.8 


NO 


< 14.0 


{Ee) IH 


12.6-16.2 


— 


14.5-15.2 


NO 


— 


14.4-15.2 


(^e) NH 


— 


14.9-15.2 


14.9-15.2 


16.4-16.8 


14.7-15.2 


14.8-15.2 


{E.) IH 


17.9-18.1 


— 


17.8-18.2 


16.4-16.6 


_ 


16.2-16.7 


(E,) NH 


— 


17.8-18.1 


17.7-18.1 


16.0-16.8 


17.3-17.7 


16.3-16.7 


eiH 


0.7-0.8 


— 


0.6-0.8 


1.3-1.7 


— 


> 1.3 


^ NH 


— 


0.7-0.8 


0.7-0.8 


0.7-0.9 


1.4-1.6 


1.3-1.7 


Eb/D"^ IH 


2.9-3.1 


- 


2.9-3.1 


2.8-3.2 


- 


2.8-3.2 


^f,/L>2 NH 




2.9-3.1 


2.9-3.1 


2.9-3.2 


2.9-3.1 


2.8-3.2 


sin2 6*13 IH 


> 2 • lO"'' 




< 1 • 10"^ 


> 1 • 10-4 




< 2 • 10"^ 


sin^ (9i3 NH 




> 2 • 10-4 


< 2 • 10"^ 


> 1 • 10-^ 


> 8 • 10-^ 


< 3 • 10-^ 


Hier. det. 


yes 


yes 


no 


no 


yes 


no 


Hier. excl. 


> 5a 


> 5a 






> 5a 





TABLE V: The allowed parameter space at 3a CL (2 dof ) for the points 2 and 3 of Table [J by a 540 
kt WaterC detector using all channels. The symbol '-' is used when there is no allowed area. In this 
case, we have written the confidence level with which a given hierarchy can be excluded. The symbol 
'NO' is used when there are no restrictions on the parameter space given in Eq. (116p . The energies 
are in MeV and Eh/D"^ in units of 10^^ ergs/kpc^. 

hierarchy has already been established. 

In Tables |Vl |VT] and IVIII we show the allowed parameter space for the input point 2 of 
Table [H For point 2, we have chosen a smaller value of ^, making the original and z/g fluxes 
twice as large as the original flux. There is not much difference between the results of point 
1 and 2, except that the determination of (Ee) and (Ee) are better for point 2. The accuracy of 
the determination of {E^) has little dependence on the value of ^, although it increases slowly 
for larger values of ^. 

In Tables |Vl I VI I and IVIII we show the results for point 3 of Table [H For this point a weaker 
hierarchy between (E^) and (Ee) is considered, being about 10% for point 3 and about 20% 
for point 1. The main difference between these points, is the fact that in the case iO, the 
hierarchy can no longer be recognized. As mentioned above, the establishment of the hierarchy 
dependents on whether a value of (Eg) equal to (Ex) is included in the scan over parameters. 
Had we enlarged the SN parameter space in Eq. (fT6l) to include (Eg) = 18 MeV we would 
find a very small allowed area for the WaterC detector for the normal hierarchy around the 
set of values: {Eg) = 18 MeV, {E^) ~ 10 MeV, {E^) ~ 18 MeV, ^ ~ 0.8, sin^ ^13 - 10-^ and 
Eb/D"^ ~ 3.1 X 10^^ ergs/kpc^. It is not difficult to realize that this set of parameters produces 
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17.3-18.6 
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< 1.5 


> 0.7 
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NO 


^ NH 


— 


0.6-1.1 


< 1.3 


0.7-1.3 


> 1.0 


> 0.7 
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2.8-3.2 


- 


2.8-3.2 


2.8-3.2 


- 


2.8-3.1 


^f,/L>2 NH 




2.7-3.3 


2.8-3.3 
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2.7-3.3 


2.8-3.2 


sin2 6*13 IH 


> 6 • 10-5 




< 5 • 10-5 


> 6 • 10-5 




< 2 • 10-5 


sin^ (9i3 NH 




> 6 • 10-5 


< 2 • 10-5 
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< 4 • 10-5 
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no 
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no 


Hier. excl. 
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3(7 






> 5(7 





TABLE VI: Same as Table |V] but for a 100 kt LAr detector. 



total number of events in each channel only slightly different from the iO case. Improving the 
energy resolution of the detector does not help much. For point 1 the normal hierarchy is 
excluded at 95% CL and for point 3 by less than la for the WaterC detector. Nevertheless, 
as the allowed region for the normal hierarchy in this case is very small and occurring for the 
z/g and z/j, fluxes having equal temperatures, one would naturally have a strong hint that the 
true hierarchy is indeed the inverted hierarchy. An analog discussion could be performed for 
the case nO. However, as this would require all three average energies to be almost equal, we 
will refrain from this discussion, as it is physically very improbable. 

Next, we would like to discuss how the detection of Earth matter effects and shock wave 
effects can help to pin down the hierarchy in the iO scenario when allowing the Ue and fluxes 
to have identical temperatures. Let us shortly review the facts about Earth matter effects for 
SN neutrinos 23|, |2J]. For neutrinos that traverse the Earth mantle (and core) a modulation 



with known frequencies of the neutrino energy spectra may occur, depending on the hierarchy 
and the value of ^13 . If the hierarchy is normal (inverted) , the Earth matter effect for neutrinos 
(anti- neutrinos) depend on the value of Ph- In both cases, the strength of the modulation 
will be proportional to the difference in the original z/^; and (ug) fluxes and the value of Ph- 
Therefore, in the case of normal hierarchy and large values of ^13 {Ph — 0)) there should be 
a modulation of the anti-neutrino energy-spectra unless (Eg) = (Ex) and no modulation of 
the neutrino spectra. Similarly, in the case iO there should be a modulation of the neutrino 
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TABLE VII: Same as Table IVl but for a 50 kt Scintillation detector. 



energy-spectra unless (Ee) = (Ex) and there should be no modulation of the anti-neutrino 
spectra. This can clearly help to distinguish the hierarchy in case iO for point 1, as we do 
have a hierarchy between the neutrino and average energies. With detectors shielded and 
unshielded by Earth, telling us that there are Earth matter effects in the neutrino channel and 
not in the anti-neutrino channel, the case iO is clearly established. The expectation for point 
1, is a maximum difference of the neutrino flux arriving at Earth and the flux after traversing 
part of the Earth matter of about 20%, occurring for a neutrino energy of roughly 60 MeV. 
Obviously for point 3 it will be much more difficult to observe and determine the Earth matter 
effects, as the weak hierarchy between (Eg) and (E^) will make the overall strength of the Earth 
matter effects smaller. Furthermore, shock-wave effects 30|] may also help to identify the true 
hierarchy. These effects can be seen in the adiabatic region of the high (H) resonance, i.e. for 
large values of 6*13. In the case of inverted hierarchy a dip in the value of (Eg) as a function 
of time is expected, whereas for normal hierarchy the dip is expected for (Ee). Therefore, an 
observation of shock wave effects in the anti-neutrino channel will point toward the inverted 
hierarchy. Once again the amplitude of this effects decrease as the u^/ Ve and temperatures 
becomes closer, making it difficult to pin down the hierarchy for point 1. In conclusion, the 
hierarchy for case iO is likely to be established with the detection of the SN neutrinos at Earth 
alone. The confidence level with which this can be done increase as the hierarchy between the 
z/g and temperatures increase. With the complementary information on Earth matter effects 



34 



and shock wave effects it is very likely that the hierarchy can be undoubtedly established for 
point 1. 

Summarizing, the values of some or even all of the SN parameters as well as the unknown 
neutrino ones, might be determined simultaneously, in most cases. The values of C,, (Eg) and 
(Eg) are however difficult to determine with a high precision. Overall, there is not much dif- 
ference in the performance of all three detectors. The WaterC can access most SN parameters 
with a higher accuracy, but this is basically due to the larger mass and hence statistics. The 
LENA-type detector is performing slightly better than LAr (except for the determination of 
{Ef.)) for the cases al and iO. For the nO case they are almost equally good, with LAr deter- 
mining (Eg), {Ex) better and LENA doing a better job in determining ^. Therefore, generally 
the two detectors have similar performances, although the mass of the LAr detector is twice 
as big as that of LENA. This difference can be understood as the LAr detector has a rather 
weak sensitivity to the pure anti-electron neutrino flux. The LENA-type detector, on the other 
hand, has a good distribution between sensitivity to pure z/g, pure Ve and the total neutrino 
flux. 

V. CONCLUSIONS 

We have studied the prospects for extracting, simultaneously SN parameters and neutrino 
oscillation parameters from the measurements of neutrinos from the cooling phase of a galactic 
supernova in three different detectors: a megaton-scale water Cherenkov, a 100 kton liquid 
Argon and a 50 kton scintillation detector. 

In our analysis we have varied a total of seven parameters, five SN parameters: the average 
energies (-Ee), {E^) and {Ex)-, the ratio of the luminosities ^ and the overall normalization of 
the fluxes Eb/D"^; two neutrino oscillation parameters: the angle ^13 and the neutrino mass 
hierarchy. Since we considered perfect detectors, with 100% efficiencies, our analysis must be 
viewed as an estimation of the maximal performance of each experimental setup. We do not 
include Earth matter effects or shock-wave effects in our calculations but briefly discuss their 
possible implications. 

Our main results are summarized in Figs. [MZl IHlfTO] and [TTlfTSl for the WaterC, LAr and 
Scintillator detectors, respectively. We have found that SN parameters, as well as the unknown 
neutrino ones, can be determined simultaneously, in most cases. Comparing the three detectors, 
there is not much difference in their overall performance. However, the WaterC detector can 
access most SN parameters with higher accuracy (a few % in some cases), but this is basically 
due to its larger mass and statistics. 

All of the studied detectors have the possibility to determine the neutrino mass hierarchy if 
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sin^(^i3) ~ 10"^ and the hierarchy of the average energies is stronger than about 20%. They 
can at the same time determine some of the SN parameters quite well and put strong bounds 
on the value of ^13. The average energy of the u^^ and Ur species can be determined with an 
accuracy better than 5% in most of the parameter space suggested by SN simulations. 

Although the detection of several separable channels measuring different combinations of 
CC and NC processes is crucial for the determination of 6*13 and the hierarchy, there are cases 
where some SN parameters can be determined rather well even when only the main CC detection 
channel is available. 

APPENDIX A: Allowed regions for detectors at non-optimal performances 

In this appendix we show contours for a few cases, where we assume that some channels are 
either non-separable from other channels or not available at all. 

We investigate the case of a WaterC detector where the inverse beta decay channel cannot 
be separated from the CC on Oxygen. We have chosen this scenario, as indeed this separation 
might be difficult if gadolinium will not be added (in which case the neutron is not detectable), 
since both detection processes are almost isotropic. We show the al case in Fig{TH Comparing 
to the case of all channels being separable in FigJS] (solid lines) we see that the biggest difference 
is the determination of sin^ ^13 for the normal hierarchy allowed region. The upper bound on 
sin^ 6'i3 increases from 2 x 10~^ to 5 x 10~^. Besides this the other parameters have very similar 
restrictions. Moreover, the contours for inverted hierarchy are almost identical to the ones in 
Figl5] (solid lines). This is expected as for the inverted hierarchy the inverse beta decay channel 
dominates and therefore whether or not it can be separated from the CC-0 does not have a 
large effect. This is also the reason that the scenario iO does not change much if the IB and 
CC-0 cannot be separated. The only visible difference is a slightly worse determination of (Ef.), 
but this variable can not be determined very well in either scenarios. For the scenario nO, the 
lower bound on sin^ ^13 jumps from 1 x 10^"^ to 6 x 10"^. Overall, the separation of the IB and 
CC-0 channel have little impact on the detector performance, besides a less restrictive bound 
on the CHOOZ angle. 

Next we study the LAr detector in the case that the neutral current on Argon is not available. 
This NC channel is not accompanied by the detection of the positron or electron and might 
therefore be difficult to trigger, especially if the detector is not well shielded. We show the case 
where the input scenario is nO in FiglTSl Comparing FigHS] to FigHO] (solid lines) we see that 
again the biggest impact is a worsening of the determination of sin^ 61^. In this case the impact 
is severe, lowering the upper bound from 10""^ to 7 x 10"^. Other parameters are only mildly 
affected. For the case iO, the normal hierarchy cannot be ruled out with the 3 channel scenario. 
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sln'(0,3) 




2 2.5 3 3.5 4 

E,/D' X 10" 



540 kt Water C 

Any hierarchy assumed (a1) 

— IH 3 (T (inseparable IB and CC-0) 

— NH 3 (T (inseparable IB and CC-0) 



FIG. 14: Sensitivity of a 540 kton WaterC detector assuming true SN parameters as in point 1 of 
Table HI for any hierarchy and sin2(6'i3) = 10-6 (Pff ~ 1). We show 3cr CL contours (2 dof) usine all 
4 detection processes (IB + CC-0 + ELAS+ NC-0), but assuming that the IB and CC-0 cannot be 
distinguished. The contours for the normal hierarchy (NH) is marked by the dark (red) and for the 
inverted hierarchy (IH) is marked by the light (cyan) solid line. 
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FIG. 15: Sensitivity of a 100 kton LAr detector assuming true SN parameters as in point 1 of Table 
m for normal hierarchy and sin^(^i3) = 10"^ {Ph — 0). We show 3a CL contours (2 dof) using 3 
channels (feCC+ v'eCC + ELAS), assuming that the NC on Argon is not available. The contours for 
normal hierarchy is marked by the dark (red) and for inverted hierarchy is marked by the light (cyan) 
solid line. The inverted hierarchy is ruled out by more than 5cr. 
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sln'(0,3) 




2 2.5 3 3.5 4 

E,/D' X 10" 



50 kt Scintillator 

Normal hierarchy assumed (nO) 

— IH 3 (T (inseparable CC on Carbon) 

— NH 3 (T (inseparable CC on Corbon) 



FIG. 16: Sensitivity of a 50 kton Scintillation detector assuming true SN parameters as in point 1 of 
Table m for normal hierarchy and sin^(0i3) = 10"'^ {Ph — 0). We show 3a CL contours (2 dof) using 
all 6 detection processes (IB + feCC+ UeCC + u — p + NC + ELAS) but assuming that the fgCC 
and i^eCC cannot be distinguished. The contours for normal hierarchy is marked by the dark (red) 
and for inverted hierarchy is marked by the light (cyan) solid line. The inverted hierarchy is ruled out 
by more than 5ct. 
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as a small 'fake' region around {E^) = 17 MeV is allowed for small values of sin^6'i3. This 
also leaves sin^ ^13 unrestricted for the inverted hierarchy, unless one constrains (Eg) to be less 
than about 16.5 MeV, in which case the restriction is roughly as the case with 4 channels. For 
the input scenario al, again the biggest impact is a drop in the upper bound on sin^ 6*13 from 
2 X 10"^ to 6 X 10~^. We conclude that the NC channel on Argon is an important factor for 
determining the neutrino parameters: the hierarchy and in particular the value of the CHOOZ 
angle. 

Finally we investigate the Scintillator detector assuming that the neutrino and anti-neutrino 
CC reactions on Carbon cannot be separated. We show the case of the input scenario nO 
in FigJT6l Comparing Fig{T6] to Fig{T3] (solid lines) we see that there is hardly any visible 
difference. This is also valid for the two other cases (al and iO). The reason for this is that the 
dominant inverse beta decay channel is still a separate channel and thus gives good restrictions 
on the Ue flux, which in turn makes a de facto separation of the Ug-CC and Ue-CC channels 
when doing the fitting. We conclude that the separation of the CC channels on Carbon only 
have a minor impact on the performance of the Scintillator detector. 
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